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摘 E 本 文 将 经 典 的 形状 灵敏 度 分 析 方 法 与 一 种 改进 的 水 平 集 方法 相 结合 ， 给 出 了 Navier-Stokes 问题 
形状 优化 的 一 种 新 方法 。 该 算法 是 在 固定 的 Euler 网 格 .上 进行 计算 且 在 优化 过 程 中 不 需要 对 水 平 
集 函 数 进 行 重新 初始 化 ， 从 而 可 以 有 效 地 节省 计算 时 间 。 数 值 算 例 说 明 该 算法 是 稳定 、 高 效 的 。 
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作为 最 优 控制 的 一 个 分 支 ， 流 体力 学 中 的 形状 最 优化 问题 很 早 就 得 到 科技 工作 者 和 工程 师 
们 的 关注 ， 并 已 取得 许多 非常 好 的 成 果 [" 引 。 在 形状 最 优化 的 问题 中 ， 最 常用 的 工具 是 灵敏 度 
分 析 。 但 是 在 通常 情况 下 ， 基 于 灵敏 度 分 析 的 优化 算法 在 计算 过 程 中 都 需要 网 格 的 重新 划分 ， 
这 使 得 计算 量 非 常 大 。 

FH Osher 和 Sethian 提出 的 水 平 集 方法 外 ， 最 近 被 引入 到 形状 最 优化 问题 的 研究 中 ， 并 取得 
了 很 多 非常 好 的 应 用 巨 引 。 与 传统 的 形状 最 优化 方法 不 同 ， 水 平 集 方法 是 在 固定 的 欧 拉 网 格 上 
进行 的 ， 即 在 计算 过 程 中 网 格 只 需 剖 分 一 次 ， 从 而 能 有 效 的 节省 计算 时 间 。 

在 水 平 集 方法 的 实施 过 程 中 ， 移 动 的 波 前 或 界面 是 用 水 平 集 函数 的 零 水 平 集 来 表示 的 。 这 
要 求 水 平 集 函 数 在 演化 的 过 程 中 要 保持 为 有 向 距离 函数 ， 但 这 一 性 质 是 很 难保 持 的。 一 种 常用 
的 修正 方法 是 在 演化 的 过 程 中 对 水 平 集 函数 进行 不 断 的 重新 初始 化 。 但 重新 初始 化 的 过 程 是 非 
常 复杂 也 非常 费时 的 ， 而 且 对 边界 有 很 强 的 依 束 性。 最近 Chan 等 学 者 提出 了 一 种 基于 变 分 原 
理 的 水 平 集 方法 ， 这 一 方法 在 水 平 集 函 数 演化 的 过 程 中 不 再 需要 对 水 平 集 函 数 进行 重新 初始 
化 Bq。 从 而 可 以 进一步 提高 计算 效率 。 

本 文 在 前 人 工作 的 基础 上 ， 特 别 是 水 平 集 方法 在 图 像 处理 及 结构 力学 形状 优化 方面 的 应 
用 ， 给 出 了 Navier-Stokes 问题 形状 最 优化 的 一 种 新 算法 。 在 这 个 方法 中 ， 我 们 将 传统 的 形状 
灵敏 度 分 析 与 改进 的 水 平 集 方法 相 结合 ， 该 算法 在 优化 过 程 中 不 需要 网 格 的 重新 剂 分 ， 也 不 需 
要 对 水 平 集 函数 进行 重新 初始 化 ， 从 而 有 效 地 减少 了 运算 工作 量 和 计算 时 间 。 后 面 提供 的 数值 
算 例 也 验证 了 这 一 点 。 
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2 ”优化 问题 模型 及 形状 灵敏 度 分 析 


BN e R? 是 一 个 具有 光滑 边界 T := 80 = T。UT4 的 有 界 开 区 域 ， 考 虑 目的 二 维 不 可 斥 
缩 定 常 Navier-Stokes 问题 


~vAut+(u-V)u+Vp=f, 在 9 A, 


divu = 0, 在 9 A, (1) 
u=0, 在 re E, 
u=g, 在 Ta 上 ， 


其 中 v = 1/Re (Re 为 雷诺 数 ) 表示 动力 学 粘性 系数 ，u 表 示 速 度 ，p 表示 压力 ，f 为 体积 力 。 

设 速度 ua 的 检验 函数 空间 为 Hi (Q) = {v € HQ), vlr. = 0}。 压 力 项 p 的 检验 函数 空间 
为 Zi(Q) = {p € L?(Q), fo pda = 0}。 为 了 定义 Navier-Stokes 方程 的 弱 形 式 ， 引 入 下 面 的 双 线 
性 形式 


a(u, v) = v(Vu, Vv) = vf Vu: Vvdr, Vu, v€ H'(Q) (2) 
Q 
和 
b(v,q) = -f qdivvdz, Yv e€ H'(Q), (3) 
Q 
以 及 三 线性 形式 
c(w; u, v )= fw Vu- vdz = > Le wj =~ )vidz, Vu, v, weHI(D). (4) 


由 Green 公式 ， 可 以 得 到 下 面 的 弱 形式 ;, SK (u, p) € HG(Q) x LQ), 1G 
a(u, v) + c(u; u, v)+Tbvp)=(fv)，VveHID)， 
b(u, g) = 0, Vq € L3(Q), (5) 
(u, Ss)ra = (g, S)ra, vse H-1/2(Ta), 


其 中 (yp, RR HTT.) := (HTA AHT.) LEIBH. ER, H Sobolev 空间 的 
KA SE FELL Holder 不 等 式 ， 我 们 可 以 得 到 :a(,:) MoC -) PWE HN) x H (QAH (N) x 
L2(Q) 上 连续 ， 而 c(,…) 在 H1(Q)xH1(0)xHi(0) 上 连续 。 从 而 可 知 Navier-Stokes 方 程 (5) 的 
解 是 存在 且 唯 一 的 。 
目标 泛 函 的 一 般 形 式 为 
(9) = | itu vadr (6) 


形状 最 优化 问题 的 一 般 提 法 为 : 求 上 在 满足 Navier-Stokes 问 题 (5) 的 情况 下 使 得 目标 泛 
PAI (6) 达到 最 小 。 

为 了 应 用 梯度 算法 ， 我 们 将 使 用 一 些 经 典 的 形状 灵敏 度 方法 的 结论 。 设 有 界 开 集 Q CD C 
R2 为 一 参考 区 域 ， 其 中 卫 为 包含 所 有 可 行 区 域 的 于 作 区 域 ， 考 虑 如 下 形式 的 区 域 


Rb) = (Id + h)(Q) := {z + h(x), z € Q}. 


其 中 h e WbO(R?, R?) 为 一 个 小 的 摄 动量 ，1d ASAT, AQ) 是 对 初始 区 域 Q 的 一 
个 摄 动 ， 而 (1d 十 hb) 是 R2 上 的 一 个 微分 同 肛 。 
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目标 泛 函 了 (Q) 在 QQ 的 形状 导数 定义 为 WH%(R?,R?) 上 在 0 点 处 的 Fréchet 导数 为 
J ((Id + h)(2)) = 7(0) + J'(M)(h) + o(h), (7) 


其 中 rr(Q) 在 机 bce(R2,R2) 上 是 线性 连续 的 。 
由 Hadamard 结构 定理 可 知 ， 方 向 导数 .了 '(Q)(h) 只 依赖 于 边界 下 上 的 法 向 迹 h.n， 即 ， 
引 理 10 设 Q 为 具有 光滑 边界 的 开 区 域 且 ij(u Vu) € W11(R2)， 则 目标 泛 函 (6) 在 Q 上 
可 导 ， 且 其 导数 为 


J'(2)(h) = [intos Vu) )dz = [ro - n(x)j(u, Vu)ds, (8) 


对 任意 的 h e W!°(R?,R?). HEP n ARRANA E. 
本 文 研 究 关 于 Navier-Stokes 问题 (5) Mig) 42517) lal, BU ERZE 


min F(9) = L | Vulde, (9) 


其 中 也 为 包含 所 有 容许 区 域 Q KR, u 为 方程 组 (5) 的 解 。 

由 Green 公式 及 引 理 1， 可 以 得 到 。 

定理 1 ” 设 Q 为 具有 光滑 边界 的 有 界 开 区 域 ，h e W)°(R?,R?). fant 及 (5) 的 解 吕 足 
够 光滑 ， 如 设 f e H9), ue H?N), WEERA (9) 的 形状 导数 为 


aonr S (EEY - SE Fe Jas, (10) 


其 中 (u, p) 为 方程 组 (5) 的 解 ， 而 (w, q) Ei FRIN EARE 
—vAw — (uVw + (Vw)u) + Vq = —2vAu, #24, 
divw = 0, 在 Q 内 ， (11) 
w=0, 在 工 E. 


3 ”改进 的 水 平 集 方法 
D C R2 为 包含 所 有 可 能 区 域 的 9 容许 形状 区 域 ， 设 存在 一 个 函数 ol) BK ERR, 

满足 

gz) >0, VreEQ, 

g(x) =0, Vee Ti=0Q, 

d(x) <0, Vee D\Q, 
单位 法 向 量 以 及 表面 曲率 上 < 分别 为 

V9 


= Ve -yonadir ne 
= Wa 和 k=V n= div( =): 


在 优化 的 过 程 中 ， 区 域 的 边界 由 水 平 集 函 数 g(z) 的 零 等 值 面 来 表示 。 水 平 集 函 数 b(z) 可 以 
通过 求解 如 下 的 Hamilton-Jacobi 方 程 得 到 


n 


sy - VIV¢(2,t)| =0. (12) 


如 上 所 述 ， 在 优化 的 过 程 中 ， 要 保持 水 平 集 函 数 为 一 个 距离 函数 。 我 们 知道 距离 函数 必需 
满足 性 质 |Vew(z)| = 1; 反之 ， 如 果 一 个 函数 9(z) WE |VE(x)| = 1， 则 p(x) 在 相差 一 个 常数 的 
情况 下 也 是 一 个 距离 函数 。 所 以 我 们 引入 如 下 的 积分 式 来 衡量 一 个 函数 8(z) 是 否 为 ER A 
的 距离 函数 


1 
Doa) = | 5(IV9l -1)°ae. (13) 
D 
Be O(a) 在 区 域 卫 的 边界 上 的 法 向 导数 为 零 ， 即 
ago op 
“m TO 在 OD L, (14) 


则 全 (8(z)) KT O(a) E y HER Gateau 导数 为 
ID(9(z)) 四 
(ED y) =- f (ae - eae, 
也 定理 1， 形 状 导 数 可 以 写 为 如 下 的 形式 
aJ(Q)h = | Vh- nds. 
J (0) f nds 


因为 在 形状 导数 (10) F, uw UR n eee KMD LAM, PAT ERS DHA 
仪 是 边界 工 上 定义 下 降 方向 h = ~Vn。h 的 法 向 分 量 h.n = —V 可 以 做 为 Hamilton-Jacobi 方 
程 (12) 的 对 流 方 向 。 在 我 们 的 算法 中 ， 水 平 集 函数 由 下 面 改进 了 的 Hamilton-Jacobi 方 程 求 得 


PAE.) _ vivg(e,t)| - MAe(2,t) — r) =0, (15) 


其 中 入 为 平衡 (A9(z,t) 一 K) 项 影响 的 一 个 正常 数 。 在 方程 式 (15) 的 左 端 ， 第 二 项 使 得 水 平 集 
函数 的 零 等 值 线 趋向 于 区 域 的 边界 ， 而 第 三 项 使 得 水 平 集 函 数 p(x,) 在 优化 过 程 中 保持 为 一 个 
距离 函数 。 所 以 在 整个 过 程 中 重新 初始 化 不 再 是 必要 的 了 。 


4 ”数值 算法 及 算 例 


在 上 面 理论 分 析 的 基础 上 ， 可 以 得 到 本 文 所 采用 的 算法 。 
1) 给 定 初 始 区域 0o， 并 相应 地 初始 化 水 平 集 函数 如 (z); 
2) MTK>OMTAR, HAWK: 
2.1) 由 方程 组 (5) Al (11) 分 别 求 得 状态 变量 u AGERE we ait (10) 式 求 得 目 
标 谤 函 的 形状 导数 ; 
2.2) 求解 改进 了 的 Hamilton-Jacobi 方 程 (15)， 得 到 新 的 水 平 集 函 数 +1。 新 的 区 域 边 
RV eta = OQ 由 水 平 集 函 数 orti 的 零 等 值 面 确定 。 
在 上 面 所 提 算 法 中 ， 不 需要 考虑 水 平 集 函 数 的 重新 初始 化 过 程 。 
考虑 经 典 的 流 线 体 的 优化 问题 。 该 算 例 曾 经 在 文献 [1,3] 中 被 研究 过 。 设 计 区 域 如 图 1 所 
未 ， 流 体 以 恒定 的 速度 从 左 向 右 ， 从 靠近 区 域 边界 的 部 分 流 过 ， 其 余部 分 为 需要 优化 的 部 分 。 
容许 形状 集 为 
D={9: QCD, |Q] = 7P] }, (16) 
其 中 | :| 表示 在 积 ，7 为 面积 比 ， 是 0 到 1 之 间 的 一 个 已 知 常数 。 为 了 与 文献 [1,3] 中 的 计算 结果 
做 比较 ， 我 们 取 雷 诺 数 Re = 10。 入 取 值 都 为 1。 优 化 过 程 中 都 没有 对 水 平 集 函 数 进行 重新 初始 
化 。 
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Al: 流 线 体 的 设计 区 域 


作为 比较 ， 这 里 考虑 了 两 种 情形 ， 一 种 是 体积 比 为 7 = 0.75， 一 种 是 体积 比 为 7 = 0.9。 
图 2 至 图 5 是 所 得 到 的 计算 结果 ， 其 中 左边 一 栏 是 体积 比 为 7 = 0.75 的 情形 ， 右 边 一 栏 是 体积 
比 为 7 = 0.9 的 情形。 图 2， 图 3 分 别 给 出 了 得 到 最 优 形状 地 对 应 的 流 场 图 。 图 4， 图 5 是 得 到 
最 优 形 状 时 水 平 集 函 数 的 表面 视图 。 图 6 为 目标 泛 函 在 每 一 步 迭 代 时 的 值 。 从 计算 结果 可 以 看 
出 ， 用 本 文 所 提 算 法 得 到 的 结果 与 文献 1,3] 中 所 得 的 结果 是 一 致 的 。 


t 


ST a a ee op} - 


Oth — -> 


Mennt te 


o 02 0.4 0.8 0.8 t 


图 3: ”最 优 形状 对 应 的 流 场 (y=0.9) 


图 4: KPRM RERA (7 = 0.75) Ad: 水平 集 函数 四 的 表面 视图 (7 = 0.9) 
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5 


图 6: 流 线 体 优化 问题 中 目标 泛 函 的 收敛 过 程 


结论 


本 文 提出 了 一 种 Navier-Stokes 问题 形 状 优化 的 新 算法 。 与 以 前 的 算法 相 比 ， 本 算法 存在 如 


下 优点 : 算法 容易 实施 ， 并 且 是 在 固定 的 Euler 网 格 上 进行 计算 ， 也 不 需要 重新 初始 化 ， 因 而 
运算 量 大 大 减少 ; 对 初始 形状 和 拓扑 结构 依赖 性 不 是 特别 强 ， 这 使 得 在 实际 计算 中 初始 形状 的 
选取 和 水 平 集 函 数 的 初始 化 更 容易 。 数 值 实验 结果 也 表明 了 算法 的 稳定 性 、 高 效 性 和 对 初始 形 


状 的 不 敏感 性 。 
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An Improved Level Set Method for the Shape Optimization of 
Navier-Stokes Problems 


DUAN Xian-bao, QIN Xin-qiang 


(School of Sciences, Xi’an University of Technology, Xi’an 710048) 


Abstract: Based on the classical shape sensitivity and the variational principle, a new level set 
method is proposed for the optimal shape control of the steady-state Navier-Stokes fluid flow. The cost 
of this method is moderate since the shape is captured on a fixed Eulerian mesh. Furthermore, unlike 
the classical level set method, the re-initialization procedure is not necessary during the optimization 
process. Promising features of the proposed method are illustrated by several numerical examples. 
Keywords: shape optimization; sensitivity analysis; Navier-Stokes problem; level set methods 
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